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Abstract 

In the framework of the Schwinger Model for percolating strings we establish a 
general relation between multiplicity and transverse momentum square distri- 
butions in hadron-hadron and heavy ion collisions. Some of our results agree 
with the Colour Glass Condensate model. 



Experimental data from RHIC (Relativistic Heavy Ion Collider) show very interesting 
features concerning particle rapidity densities and transverse momentum, p?, distributions 
[TJI2]. They exclude high particle densities, expected in naive multicollision models [3], as 
well as fast growing values for (px) as a function of energy, expected in naive perturbative 
QCD models [I]. Physics seems to remain classical and, essentially, non-perturbative. 

Multiparticle production is frequently described as resulting from multiple collisions at 
the parton level and, in the case of nucleus-nucleus collisions, also at nucleon level, with 
formation of colour strings stretched between the projectile and the target, which decay 
into other strings that subsequently hadronize into the observed hadrons j^j- There are 
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long strings in rapidity, valence strings, associated to valence quark (diquark) interactions, 
and short strings in rapidity, centrally produced (sea strings) associated to interactions 
of sea partons, mostly gluons. In a symmetrical AA collisions, with Na participants from 
each nucleus, the number of valence strings equals the number of participants, as in the 
wounded nucleon model [6 , while the number of sea strings behaves roug hly as N s » N% 3 
[Zj, increasing with the energy 

We shall adopt for the mechanism of particle production the Schwinger model mecha- 
nism as developed in [HI El- In particular, the particle density and transverse momentum 
square will be considered proportional to the field (and the charge) carried by the string. 

In multicollision models, many strings are produced, the number increasing with en- 
ergy, atomic mass and centrality. If the strings are identical and independent, and ap- 
proximately align with the collision axis, we have, for the rapidity particle density, dn/dy, 
and for the average of the square of the transverse momentum, (py), 

Ty = (1) 

<Pt) = A, (2) 

where N s is the number of strings, fi\ is the single string particle density and p\ the 
average transverse momentum squared of the single string. 

If the strings fuse in a rope [0], the colour randomly grows as \/W s and we have 

dn 1 , T 

-.N s n u (3) 



dy y/N s 

(P 2 T ) = P? V 7 ^- (4) 

In the situation of a hadron-hadron or nucleus-nucleus central collision, the strings 
overlap in the impact parameter plane and the problem becomes similar to a 2-dimensional 
continuum percolation problem |10j . If the strings are randomly distributed in the impact 
parameter plane then, in the thermodynamical approximation ^T] , the overlapping colour 
reducing factor is given by 

ll - e -v 

F(v) = (5) 
where rj is the transverse density percolation parameter, 

H^) 2 ;v, < 6 > 

where nr 2 s is the string transverse area and 7rR 2 the interaction transverse area. We thus 
have 

dn 

— = F^Nm, (7) 
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Equations similar to (jZJ and (jHJ) were written in ^Tj. As with rj — > (low density limit) 
F{rj) — > and with 77 —>■ 00 (high density limit) -F(^) — > 1/ y/r), the behaviour of relations 
dU and (0), and © and (@j) is recovered from © and ©• 

We shall now discuss the consequences of (JZ|) and (JHJ). Two straightforward results 
follow: 

i) slow increase of particle density with energy and saturation of the normalised particle 
densities as N s increases 

As the number of strings, N s , increases with energy, at large energy 77 also increases 
and 

F(V) ~ (9) 



Nl /2 ni . (10) 



which means, (jZJ), 

(in /i?" 

Instead of growing with iV s , as one should have naively expected with independent strings, 
(HJ), the density grows more slowly, as N^ 2 . 
On the other hand, as 

N s ^Nf , R^R.N 1 /, (11) 
where Ri is a quantity of the order of the nucleon radius, 

1 dn f Ri\ / \ 

— — ^ [ — )ni 12 

N A dy \r s J v ' 

tends to saturate as Na increase. Both behaviours (JZJ) and ()12)1 were confirmed by data 

m 

The saturation, in our framework, is a consequence of string percolation [I2j- At the 
level of QCD it can be seen as resulting from low- a; parton saturation in the colliding 
nuclei 



ii) a universal relation between dn/dy and (pr) 
For large density, Eqs. ((Zj) and (jHJ) become 



(Pi) = (£) A'," 2 ??, (14) 



and, eliminating iV, 



with 



1/2 



( ^ ) = C ^ay (15) 



, — x 1/2 

-.2 X ' 



fH • as) 



i?l / \ 77-1 
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A relation of this type, 



1 dn 



was obtained, in the framework of the Colour Glass Condensate (CGC) model |14j . in 
[TH] . Our formula (fT4*|) includes not only the functional dependence, but, as well, the 
proportionality factor c. 

We can make an order of magnitude estimate of the proportionality factor c. In the 
dual string model r s pa 0.2 fm [TU1 I16j. Ri should be of the order of the proton radius (~ 
1 fm) and for the string charged particle production parameters one has pi ~ 0.3 and fii pa 
0.7, as observed from low energy data and (p\/ni) l l 2 ~ 0.35. The proportionality 
factor is then pa 0.07 to be compared with 0.0348 for pions and 0.100 for kaons |TK] . 



In the comparison with data we shall identify \j (p T ) with (pr) and y pi with p\ (this 
overestimates the average values of (px) and pi). 

We have just considered the high n limit. In the low density end, which means low 
energy and peripheral collisions, we have just valence strings and (px) — > pi ~ 0.3 GeV. 
This is, in practice, the value of (py) in pp collisions at low (y/s < 10 GeV) energies. 

By putting these two limits together, we arrive at the formula obtained in |15j . but 
now with all the parameters theoretically constrained: 



_ . r s 1 1 dn 

In Fig. ^ we compare Eq. (fT%|) with data. The agreement is not perfect, but there is an 
indication that some truth exists in CGC and string percolation models. 

In the next step we make an attempt to generalise our results and to relate the (nor- 
malised) transverse momentum distribution f(p T ) to the multiplicity distribution P(n), 
in hadron-hadron and nucleus-nucleus collisions. 

We work in the large r\ limit and start by changing the notation, and write 

N = aNl /2 , (19) 

with 

a = R/r„ (20) 

such that iV has the meaning of the number of effective strings (mostly sea strings or 
ropes). If n particles are produced 

n = Nm, (21) 

see also Eq. (|13|). This effective number N takes into account percolation effects in the 
sum of colours of the N s individual strings. 

Let P{N) be the probability of producing N effective identical strings and p{rii) the 
probability of producing particles from the z-th string. We then have 

P(n) = I P{N) Y[p(ni)dni8 ( n - ^ nA dN (22) 
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Figure 1: (pr) vs multiplicity density in pp collisions (where Na = 1) at 1800 GeV ^H] (open circles) 
and in central Au+Au collisions at 200 AGeV 19 (filled squares). Solid lines represent Eq. 1)18(1 with pi 

adjusted separately to each species. 



In (12 2|) . as the colour percolation effects were absorbed in N, we treated the effective 
strings as independent (see [2T)|). 

Regarding transverse momentum distributions, the natural generalisation for (fT4*j) is 
to write 

N 



P 2 t = —,vl (23) 



and for the distribution itself 



F(Pt) = J P(N)f(pl)5 U T - ^ dpjdN. (24) 

In this case, for a given F(p^) contribute all effective strings with f(pf), such that p\ 
satisfies (|2~3|) . As all strings are assumed equal, f{p\) is representative of any string. 

In order to construct P(n) and F(p^), one of course needs the elementary string 
distributions p(ni) and f(pf) and the distribution P{N) of effective strings. Concerning 
the p(rii) distribution, it should be Poisson or close to Poisson type (as seen in e + e~ 
at low energy [21] )• The p\ distribution in the Schwinger model is an exponential in 
— p\. The P{N) distribution contains the nucleonic and the partonic structure of the 
colliding particles and the combinatorial factors of Glauber-Gribov calculus; its shape is 
investigated in [2~2~| . 

Our objective here is not to solve Eqs. (j22J) and (|2*4*jl . but simply to try to relate P(n) 
to F(p 2 T ). In view of that, let us proceed by calculating the (n q ) and (p^) moments of 
the distributions (J22|) and ()24|) . respectively. The calculations are straightforward, but 
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lengthy in the case of multiplicities (see, for example, [201 ) - I n this case, to simplify, we 
shall assume 

p(n 1 ) = 5(n 1 - ni). (25) 

It has been shown, sometime ago, that this approximation in hadron-hadron and nucleus- 
nucleus collisions is very reasonable |20| . 
We have then for the moments: 



(n q ) = (N q )n\. (26) 

It is clear, because of ([25)1 . that all fluctuations come from fluctuations in the number 
of effective strings. For distribution 



Eqs. ()26|) and (|2~7jl are the natural generalisation of f)13|) and (I14J) . As before, the moments 
of the effective string distribution can be eliminated by dividing ()27|) by ()26|) and a relation 
between (n q ) and (p^?) established. But one can now do better and eliminate the strongly 
model-dependent parameter a = R/r s , eq. fl2()l) . If one writes the KNO moments 



C f = j^y n > g = i,2,..., (2f 



the parameter a disappears. By using a capital C for final distributions KNO moments 
and a small c for single string distributions KNO moments, our final result can be written 

as 

^ = %- (29) 

c ? cf 

This equation, as mentioned before, is strictly correct only for = 1. 

It is not easy to check Eq. (|2*5j) accurately, as most experiments can only measure px > 

{Pt), but one can nonetheless attempt a somewhat rough comparison. In the Schwinger 

p 2 

model the pt distribution is Gaussian, which means c q T = q\. If the final p? distribution is 
also a Gaussian, then one obtains C™ = 1, which is not a good approximation. If the final 

_ ^ 2 

Pt distribution is an exponential, which is closer to reality then C q T = (2g+l)!/(3!) , 
and we obtain, for instance, = 5!/(3!) 2 2! pa 1.66. This is to be compared with the 
experimental value C 2 n pa 1.3 at ^ = 200 GeV [21]. 

Finally, the main point we want to make with (|29|) is that multiplicity and transverse 
momentum distributions are deeply related: one should remember that, in general, from 
the Cq moments one can construct the distribution in KNO form, (x)P(x/(x)). 
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